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The region of small transverse momentum in qq~ and gg— initiated processes must be studied in 
the framework of resummation to account for the large, logarithmically-enhanced contributions to 
physical observables. In this paper, we will calculate the fixed order next-to-leading order (NLO) 
perturbative total and differential cross-sections for both a Standard Model (SM) scalar Higgs 
boson and the Minimal Supersymmetric Standard Model's (MSSM) pseudoscalar Higgs boson in 
the Heavy Quark Effective Theory (HQET) where the mass of the top quark is taken to be infinite. 

■^4- ' Resummation coefficients B g 2 ^ , C g 2 g for the total cross-section resummation for the pseudoscalar 

f*"^ ■ case are given, as well as C gg for the differential cross-section. 
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I. INTRODUCTION 



f\J ■ The discovery of one or more Higgs bosons is the central research interest for high energy physics programs at hadron 
\ colliders around the world. Beyond the phenomenology of the Standard Model (SM) Higgs boson, the Minimal 
ON i Supersymmetric Standard Model (MSSM) which is a special case of the Two Higgs Doublet Model (2HDM) is of 
particular interest to theorists. For a review see Ref. 0,0. 

Very recently, a new central value for the top quark mass was reported^- This changed the exclusion limits on 
the SM Higgs boson, putting its central mass value from precision elcctroweak fits at 117 GcV/c 2 . This is exciting 
because this value is above the exclusion limits from the LEP direct searches which exclude the mass of a SM-like 
Higgs boson below approximately 114 GeV/c 2 £|. In the MSSM there are five physical Higgs bosons; a light and a 
heavy scalar (h°,H°), two charged scalars {H ± ), and a CP-odd pseudoscalar (^4°). The mass of the pseudoscalar is 
excluded^ from being lighter than 92 GeV/c 2 . The ratio between the vacuum expectation values (vevs) of the two 
neutral Higgs bosons of the MSSM is defined as tan/3 = V2/V1. For m top = 174.3 GeV/c 2 , 0.5 < tan /3 < 2.4 has been 
Oh, excluded by the LEP Higgs searches. The bounds on tan /3 will change as the central value of the top mass changes. 

A full analysis of the tan/3 exclusion bounds for the new top mass of m t0 p = 178.0 GcV/c 2 is not yet available. In 
. • , this paper, we will leave tan /3 = 1 so that the pseudoscalar results can easily be scaled by an appropriate number of 
^ ■ interest to the reader, and the mass of the Higgs bosons will be set to M$ = 120 GeV/c 2 , where $ is the Higgs boson 
of interest. 

In the context of resummation, the literature has focused on the scalar Higgs boson|rl I7I IM Ifjl Hfl ITU IT^ . H3L ITiL [TtI 

III IH 13 E E3 

This paper will provide resummation coefficients for the pseudoscalar Higgs boson for the total 
cross-section and differential distributions. Our calculations are done in the Heavy Quark Effective Theory (HQET) 
where the mass of the top quark is taken to be infinite. The role of the bottom quark in pseudoscalar production 
becomes dominant at large tan/3. In order to correctly take the bottom quark into account at this order, massive 
resummation coefficients will have to be determined. This will be reserved for another discussion. 

In Section ^ we will introduce the Heavy Quark Effective Theory (HQET) in which our calculations were per- 
formed. In Section IlIII we will introduce our resummation conventions and present our new results for pseudoscalar 
resummation. Finally, in Section llVI we will present our numeric results for the differential distributions for the SM 
scalar and MSSM pseudoscalar Higgs boson at the Large Hadron Collider (LHC) and Tevatron in the HQET. 
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II. HEAVY QUARK EFFECTIVE THEORY 



Higgs phenomenology in QCD lends itself well to the use of Heavy Quark Effective Theory (HQET). When only 
the top quark is considered in calculations, it is possible to replace the top quark loops by an effective vertex when 
the other quarks are ignored. The role of the bottom quark in Higgs physics has recently been examined in great 
detail0, but will not be included in this paper. The Lagrangian that describes this effective vertex can be 
derived from the gg — > $ (where $ is a Higgs boson of interest) trian gle diag ram |23l I24I I2H |2(| and letting the mass 
of the quark become infinitely heavy at the end of the calculation [271 l28l l29| . The Higgs of interest could be the SM 
scalar Higgs, or the pseudoscalar Higgs of the MSSM. In principle, the light and heavy scalar Higgs of the MSSM can 
be included in this formalism in place of the SM Higgs by multiplying by the appropriate coupling factor. However, 
supersymmetric corrections within the MSSM will not be included in this paper, only SM QCD corrections will be 
included. 

The HQET method allows for the inclusion of the 0(e) terms which are important for deriving resummation 
coefficients. This program leads to the effective Lagrangian in d = 4 — 2e dimensions for a scalar Higgs boson 



£f ff = -\g H HG^G a ^(^-Yv(l + c) (1) 

where gn = a s /37rw is the coupling of the effective vertex at LO, G a,)Jjl ' is the field-strength tensor for the gluons, /j, r 
is the renormalization scale, and the vacuum expectation value (vev) of the Higgs is defined as v 2 = (V^Gf)^ 1 — 
246 GeV. The effective coupling receives order by order corrections that have been calculated 0,H3 previously. The 
appearance of the top quark mass in this expression hints that the corrections to the effective coupling at higher order 
may include logarithmic corrections including the top quark mass, which is the case. Alternatively, we can define 
1/v 2 = 6422.91 pb, which is convenient in cross-section calculations. This effective Lagrangian generates effective 
vertices with two, three, and four gluons with a scalar Higgs boson. The Feynman rules for a scalar Higgs can be 
found in the literature |31j. 

When a pseudoscalar Higgs boson is considered, the effective LO Lagrangian changes due to the 75 coupling and 
can be written 

C$ = \g A oA°G^G^ (^f) V(l + e) (2) 
4 M V top / 

where g^o = a s /2irv is the coupling of the effective vertex. The G a, ^ v — 1/2 e' u /pc 'G a pa operator is the dual of the usual 
gluon field-strength tensor. It is important to note that the four gluon plus pseudoscalar Higgs vertex is absent in 
this effective Lagrangian as its Feynman rule is proportional to a completely antisymmetric combination of structure 
functions and therefore vanishes. It should be noted that this is the LO effective Lagrangian, and that a second 
operator begins to contribute at higher orders. A complete discussion can be found in Ref. |32j. The Feynman rules 
for the pseudoscalar can be found in the literature |33|. The Feynman rules listed in this reference should be used with 
5^40 = a s /27rw to avoid a spurious extra factor of 1/4. 

If we define z = M$/s where s is the partonic center of momentum energy squared, the partonic total cross-section 
for Higgs production (either scalar or pseudoscalar written here generically as $) from the fusion of two partons a 
and b can be written as a series expansion in a s as 



= <7 * (* OB M(i -*)+ (*) + ( v) A 2L*(*) + +-)> ( 3 ) 

where cr* is the LO partonic total cross-section, and the A^?^ $ (z) coefficients are higher order corrections. 

The NNLO corrections for both the scalar and pseudoscalar Higgs bosons have been calculated in the H QETlSiLIsM 
l36[ . Although the next-to- leading order corrections have been calculated in several places fA loL l33t l3ll l38l l39l l40l l4ll| . 
there are some discrepancies in the literature that we would like to clear up in this paper. For this reason, we will 
explicitly calculate the NLO corrections for the gg initial state. The importance of this particular channel, beyond 
its relevance to high-energy hadron collisions, will be addressed later. The partonic cross-section at NLO has to 
be written as the sum of the real emissions, the virtual corrections, the charge renormalization, and Altarelli-Parisi 
subtractions as follows 



<?NLO = OVeal + <7 v irt + <Jren + ^AP- 



(4) 
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Although a few very thorough treatments for the scalar case exist in the literature [33, |38(, we will re-derive them 
here to highlight differences in the pseudoscalar case, where an exhaustive treatment is missing. We will follow closely 
the discussion in these references. We will see that only the gg channel will play a role in the resummation formalism. 



A. Scalar Higgs Matrix Elements 

Suppress ing O(e) terms for now, the matrix elements for scalar Higgs production in the HQET can be written 

as 0, U Hi m OH E EH 



M(g(tf> A ) +g{v1 B ) - H) = -^-S AB (r,^^- - p^eMe v {p 2 ). (5) 

Here AB(ijlv) are the color (Lorentz) indices of the incoming gluons. When the matrix elements are squared, the 
contraction of the delta functions yields a factor of 5 ab Sab = N% — 1 = 8, the contraction of the metric yields 
V^V^ = d. This yields the squared matrix elements (before any color-spin averaging) 

1D7T V \ m top / 

In this 2 — > 1 process, it is easiest to calculate the decay width of the Higgs and convert that to a partonic total 
cross-section. We need to color and spin average the matrix elements squared noting that the gluon-gluon initial state 
must be averaged over 4(1 — e) 2 transverse polarizations. We can write the partonic cross-section in terms of the 
decay width and simply add 6(1 — z) since z = Mfj/s = 1 at threshold, so that [371 



7T 2 



a»(gg - H) = — T T(H - gg) (7) 
a s \ 2 7T /4^ 2 \ 2£ r 2 (l + e) 



7r J 576w 2 V Wtop ) 



This factor (with its e dependence) will be pulled from each of the higher order correction factors. The LO cross-section 
in the HQET starts at 0{a 2 s ). 



1. Radiative Corrections 



At the next order in perturbation theory, there are qg and qq initial state processes. However, as we will see later, 
we are only interested in the gg initial state, so we will only calculate the NLO corrections to the gg initial state 
process. 

The real contributions at NLO to the gg initial state come from the process qq — » qH which was originally calculated 
in the e — > limit in the full theory including a finite mass top quark in Ref. [42l.l43|. We can write the amplitude as, 



gipf' 11 ) + g(p B n ^ 9{-Pz a ) + h(- P6 ) (8) 

If we define the partonic (with hats) kinematic variables in terms of the Higgs momentum so that our differential 
cross-section can be written in terms of the Higgs transverse momentum as s = (pi + P2) 2 , t = (pi +P5) 2 , and 
u = (p 2 +p 5 ) 2 - The matrix elements take the symmetric form[7t WL IM I3I IM Hoi l4lj| 



\M(gg -> gH)\ 2 =-f g ° 



W r 



v- 



M% + s 4 + t 4 + u 



'top 
41 



T 2 (l + e) 



stu 



(l-2e) 



(Mfj + s 2 +i 2 + u 2 ) 2 



stu 



(9) 
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Color and spin averaging gives an additional factor of 1/256/(1 — e) 2 . By pulling out the LO cross-section from 
the expression, we can write the properly averaged matrix elements (where the overbar corresponds to color and spin 
averaging) 



\M(gg gH)\ 2 = a 



J 4 + £ 4 



",41 



stu 



(l-2e) 



(Mfj + s 2 + P + u 2 ) 2 
stu 



First, let us find the differential cross-section. The LO differential cross-section involves 2 — > 2 kinematics and can 
be written with the 4 — 2e phase space dimensions as0 



da H 1 








dt ~ 167TS 2 


K Ml J 


r(i-e) 





(10) 



As we will see in Sec. IIIII it is the small p t behavior of this expression that we will be interested in. If we insert 
the expression ut = sp 2 and drop the terms proportional to p t , we find an expression for the differential cross-section 
in the small pt limit. Once we have changed variables, we find the partonic differential cross-section in the small pt 
region (where we are suppressing the trivial rapidity dependence) 



da 



dp? 



i\i H 



(M 2 H 

(l-e)T(l-e)\p 2 



C. 



pi \ pi J 



P0 -Q 
Pt 



(11) 



We will see that this representation will make it particularly simple to extract the resummation coefficients A)^ and 
Bg 1 ^ for the differential distribution. These coefficients will be the same for the total cross-section as well. Turning 
our attention back to the expression we had for the matrix elements squared, we would like to calculate the real 
corrections for the NLO total partonic cross-section. The matrix elements for the real emission needs to be integrated 
and can be written as 13 71 



(12) 



For a 2 — > 2 process, this can be done by introducing the following parameterization for the angular integration. If 
we write the scattering angle 9 as cos 9 = 2lo — 1, this maps the 9 integration to an w integration between and 1. 
We can express the kinematic invariants as follows 



t = — J(l — z)(l — to), u = — s(l — z)uj, 
and do the phase space integrations using the following parameterization 

1 fintfYil-z) 1 - 2 * 



dPS 2 



8tt 



r(l-e) 



oj e (l — lu) e du 



(13) 



(14) 



which reduces the angular integration into repeated applications of Euler's beta function integral, where additional 
integer powers of lu and (1 — ui) are introduced from the kinematic variables t and u, 



V ' r(2 + a + (3) 



(15) 



Turning our attention back to the real emissions, we find the color and spin averaged partonic total cross-section, 
after regulating the singularity at z — 1 with a plus distribution, can be written as 



real 



xC, 



4:1th 2 . 
s 

1 1 

"2 +- 



r(i + e) 

1-2C 2 



1 + z 4 



1-z 



6(1 -z)- 
V 1 (z) 



2z 
e 



1-z 

zV (z) + +z(l - z) 



z 

2 i „4i 



1 - z + 2z z + z A V {z) 
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O(e), 



(1-zf 



(16) 
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where the plus prescription is defined as usual 



dxf(x)[g(x)} + = / dxg(x)[f(x)-f(l)}, 



(17) 



and we have introduced the common abbreviation 



VJz) 



ln n (l - ; 

1 - z 



(18) 



2. Virtual Corrections 



Next, we must calculate the virtual contributions. There are two diagrams with gluon loops (a gluon triangle 
and a four point incoming state as can be seen in Ref. [3?|l. and can be calculated directly. We can write the 
integrated virtual contribution that contributes at to the total partonic cross-section in the same fashion as the 
real emissions H S |13 S |H E E3 , 



Cvirt — CO 



T(l + e)C / 



4C 2 



6(1 -z) 



(19) 



As expected, the e 2 singularities cancel between the real emission and virtual graphs. It will turn out that the 
differences between the fixed order results for the scalar and the pscudoscalar come from different virtual corrections. 
These virtual contributions will be needed to cancel some e poles in the expression we will derive to determine the 
process dependent resummation coefficients, in particular C^g for the differential distribution. 



3. Total cross-section 



The remaining 1/e singularities must be removed to find the partonic total cross-section. The poles are cancelled in 
the charge (coupling) renormalization and the Altarclli-Parisi subtraction. Understanding the charge rcnormalization 
tells us why is was so important to have the 0(e) terms of the lowest order cross-section. We can see that the 
counter-term can be writtenl3ft|38l 



(4Z g )a , £ , ff = -^(47r) e - 1 r(l + e)A), A) = ^C A - \n f T R , 



(20) 



where = 5 since the top quark has been integrated out. These equations hold with the MS renormalization 
conditions. 

The Altarelli-Parisi subtraction factors out the soft and collinear singularities into the PDFs much like the factor- 
ization process separates the short and long distance physics in hadron-hadron scattering. This cancels the rest of 
the 1/e poles and gives us the final expression for the total cross-section 



<5nlo(<?<? —> H + X) = &(f i 5(1 — z) + I — 



^C A + 2C A ( 2 ) 6(1 -z)- Hc x (l - z f 



-2C, 



l + z 4 + {l-z) 4 



V x (z) + 2C A z 2 V (z) + (l-z) + z 2 (l - z) In 



M 

ZjJL* 



2 1 



+ 0(e). (21) 



Our expressions for the resummation coefficients show the their full color dependence. It is sufficient to notice at 
this stage that the term proportional to the 6(1 — z) in the correction can be evaluated as 11/2 + 7r 2 when we use 



C A = N C , C f 



2N„ 



(22) 
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B. Pseudoscalar Higgs Matrix Elements 



There are many reasons for our primary interest to be the pseudoscalar Higgs boson. In the MSSM, the exact 
roles the top and bottom quarks play in the differential cross-section is complicated |22j . However, in much of the 
parameter space in the MSSM, the cross-section for the pseudoscalar is larger than the lightest scalar Higgs boson in 
the MSSM. If supcrsymmetry does exist in nature, the pseudoscalar Higgs may be the first Higgs boson discovered due 
to its larger cross-section. If supcrsymmetry becomes important only at very high scales, then seeing a pseudoscalar 
Higgs would be the first evidence of supcrsymmetry in nature. This leads us to investigate in detail pseudoscalar 
resummation. 

We should also mention that because of the importance of the bottom quark in calculations involving the MSSM 
pseudoscalar (and the lightest scalar Higgs in the MSSM as well) when the parameter tan /3 is large is systematically 
ignored in these calculations. To remedy this situation, one would have to calculate the resummation coefficients in 
the full theory. In principle, these coefficients can be extracted from Ref. (3^. However, these results have not been 
published. 

The difference in the lowest order (LO) partonic cross-section of the pseudoscalar Higgs boson and the scalar Higgs 
boson can be traced to the difference in the effective couplings gu and g^o and a factor of 9/4. This can be written 
as|13 

af(gg^A°) = ^(gg^H) (23) 

= {-) 256^ ) -T^ 6 ^-^ 

Upon the expansion of the 0(e) terms, we see they do not effect the final answer as expected at the lowest order. 



1. Radiative Corrections 



The matrix elements for the production of the pseudoscalar in the HQET are slightly more complicated in d = 4— 2e 
dimensions due to the presence of the intrinsically 4-dimensional Levi-Civita tensor in the Feynman rules |33l| coming 
from the G a4iV term in the effective Lagrangian. There are several conventions for handling this problem [il . I451 litl |47| . 
We have chosen the scheme defined in Ref. 0, 0, 01 ■ 

For the radiative corrections, we can separate all the vectors into 4- and (d — 4)-dimensional components. We 
label the (d — 4)-dimensional components of the vectors with a twiddle. We can take the incoming momentum to be 
4-dimensional as a convenient choice of frame, which simplifies the results consider ably j3^. 

The (un-averaged) matrix elements can be written, 



\M(gg ^gA°)f 



a s 3 N C (N 2 - 1) f4irp 2 r 



'top 



M 



AO 



stu 



r 2 (i + e) 

2s(t 2 +u 2 ){p 3 -p 3 s-tue) 



(24) 



We can see that the e — > corrections are identical in the scalar and pseudoscalar case for the real emissions. 
We can also see that the residual difference is not only proportional to e, but rather proportional to e and the 
(d — 4)-dimcnsional component of the p 3 vector, which vanishes in the 4-dimcnsional limit. 

From this analysis, we can see that the real part of the differential cross-section for the pseudoscalar Higgs boson 
will be identical to the scalar case in Equation (|llf> with a 9/4 difference in normalization. 



2. Virtual Corrections 



The virtual corrections to the pseudoscalar have the same diagrams as the scalar at this order, but there is a slight 
difference in the result. This difference is due to the fact that the diagram with the four point gluon vertex vanishes 
due to the antisymmetry of the ggA° vertex in the effective theory. The integrated result is 



47T/^ 2 

TiJ 



r(l + e)C y 



I 



4C 2 



5(1 -z) 



(25) 
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We can see that although the pole terms are the same, the finite terms have changed slightly because of the "missing' 
diagram. 



3. Total cross-section 



When we combine our pseudoscalar results, we find that the total partonic cross-section that is identical to the 
scalar case with the exception of the small numeric difference in the 6(1 — z) term. The expression changes from 
II/6C4 + 2C4C2 — * 2(7^ + 2C A Qi in the pseudoscalar case. Here we see the factor 11/2 — > 6, which would seem to be 
a small difference, but it is mostly a coincidence of the SU3 Casimir invariants. 

The partonic cross-section can be written as 



<o(fffl - A° + X) = &f [S(l z) + 



+2C, 



l + z 4 + (l-z) 4 



\6(l-z)-^C A (l-zf 
X>i(z) + 2C A [ z 2 V (z) + (1 - z) + z 2 (l - z) In — § 



+ 0(e). (26) 



Now that we have expressions for the total partonic cross-sections for both the scalar and the pseudoscalar Higgs 
boson wc see that the only difference between the two lies in the correction proportional to 6(1 — z) at NLO and in the 
normalization. This difference in the 6(1 — z) factors is numerically small and is a s suppressed, leaving us to believe 
that the primary difference is going to be factor of 9/4 in the LO partonic cross-sections. 



III. RESUMMATION 



To introduce the machinery behind resummation |48l l49l l50l | . we need to define the hadronic cross-section. This is 
the convolution of parton distributions functions (PDFs) with the partonic cross-section 



a(S,Ml) = ^2 dxx dx 2 f a / hl (xi, fJ^f) fb/h 2 ( x 2, fJ-f) / dz z a* A ab ^(z, [i r , [i f ) (27) 

where fi r and fif are the renormalization and factorization scales respectively, and f a /t 1 is the parton distribution 
function for finding a parton a in hadron hi. We must also remember that there is a <5(1 — z) in the definition of 
the LO cross-section <r*. The minimum partonic energy fraction 2 ), m in is defined so that there is enough center 
of momentum energy to create the desired final state particles. A similar equation can be written for the differential 
distribution. 

Implicitly, the partonic cross-section contains logarithmic corrections that are formally singular at threshold (z — > 1). 
The differential cross-section contains corrections that are singular as the transverse momentum of the Higgs particle 
vanishes. They can be written in the form 



threshold - a" —, [ r Z) , recoil ~ % In 2 ™" 1 , (28) 
(1 - z) pi pi 

and various powers of these combinations. It can be seen then that the normal fixed order cross-section calculation 
diverges (in one direction or the other) at small pt due to large logarithms, and therefore it is not reliable in this 
region. The systematic way of handling these formally divergent terms at small pt is known as resummation. Because 
of this divergent behavior, one is not usually able to integrate the differential cross-section all the way down to p t = 
or to reliably understand the differential cross-section in the experimentally interesting small p t region. Resummation 
coefficients can be determined for both differential distributions and total-cross sections to address this problem. 



A. Formalism 



The resummation formalism allows the small pt cross-section to be written as a power series in both universal and 
process dependent coefficients. We write the resummed differential cross-section for a cc — ► $ process (where c in this 



8 



case represents a gluon or a quark) 

dcr rcsum 



00 b S 

db -Joibpt) f a / hl (xi,b /b) f b /h 2 (x2,b /b) 7 ^W ab (x 1 x 2 S;Q, b, 4>), (29) 



Wab(s;Q,b,, 



E 



dzi I dz2C C a{a s (bo/b),zi)Csb{a B {bo/b),z2)5(Q - z\Zis) 



do L ° 



Sc(Q,b), (30) 



where the Higgs mass Af| = Q 2 , d0 is the phase space of the system under consideration, and o)^ ^ is the lowest 
order cross-section with a cc initial state which is therefore defined at p t = 0. The constant bo is written in terms 
of the Euler-Mascheroni constant 7^ = 0.57721 . . . as 60 = 2e~ 7E . In the rcsummation formalism only the gg and 
qq initial states are needed to determine the hadronic differential distribution at small pt- The qg initial states are 
accounted for in the cross terms in the convolution. The coefficients C a b are process dependent and can be written 
as power series to be described below. Jo(bpi) is the first order Bessel function. The Sudakov form factor S c , which 
makes the integration over the Bessel function convergent, can be written as 



O 2 

A c (a s (q))ln\ + B c {a s {q)) 
1 

The coefficient functions A c , B c , and C a b can be written as power series in a s as 



{-f 


dq 2 ' 


I Jbl/b 2 


q 2 . 



n=l V / n=l V / 



(31) 



(32) 



C ab (a s , z) = 5 ab 8(l - z) + 



CO 

E 

71=1 



(33) 



The Ac , Ac 2 \ and S^ 1 coefficients have been shown to be universal. There are several conventions in the literature 
as to whether to expand in terms of a s /n or a s /2n (or even a s /47r in Ref. [3(|). We have chosen to expand in a s /n. 
It would seem that several typos exist in the literature due to this numeric expansion factor. We have derived the 
previously unknown coefficients Bg 2 \ Cgg 1 and Cg 2 ) for pseudoscalar Higgs production for the total cross-section 
and the Cgg for the differential cross-section resummation given below. Here we must stop to address a question 
of notation. It is unfortunate that we have the same notation for the resummation coefficients for both the pt 
resummation and the total cross-section resummation. It would be convenient to use a calligraphic font for the pt 
coefficients, but several authors have used this font in other contexts dealing with resummation. Therefore, we will 
put bars over the resummation coefficients for differential cross-sections even if they arc identical to the coefficients 
for the total cross-section resummation. 

To determine the coefficients for the differential cross-section, one must understand the meaning of the re- 
summation formula. We can expand Equation 129f) order by order in a s and compare to the perturbative calculation 
to read off the coefficients 01- To extract the coefficients from the perturbative results, we need to integrate the 
differential cross-section around pt = paying careful attention to the use of the Altarclli-Parisi splitting functions 
near p t = as follows 



Jo dpi 



(34) 



This expression will contain e poles and virtual corrections at the next order will be needed to be added to find a 
finite expression. This is demonstrated later in this paper. One should also be careful to use this formula with other 
quantities that show the same rapidity dependence. When we expand to 0(a s ), we can see that the NLL coefficients 
emerge as follows for a cc initial state 



Aa c5 = 1 + — 

TT 



A 



(i) 



In' 



v it 



2C). 



(35) 
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In principle, it is possible to continue this process to higher orders to obtain the needed coefficients for the differential 
cross-section. The NLO corrections to the differential cross-section are known^^^^i however the NNLO differential 
cross-section for Higgs production is currently unknown. However, the total cross-section is known to NNLO, so the 
resummation coefficients for the total cross-section can be determined to NNLL. The NLO differential cross-section 
in Ref. [l5j has been written in terms of the pt of the Higgs boson for the scalar case, and could in principle be 

- (2) 

used in part to extract the NNLO process dependent C gg coefficient for the differential cross-section for the scalar 
Higgs boson, advancing the resummed expressions ahead of the fixed order calculation. This work has not yet been 
completed. 



B. Matching 



The resummation formalism is valid in the small pt region. Fixed order perturbation theory is valid at moderate pt 
where there are no large logarithms. The process of matching allows for a smooth transition between the two regions. 
The procedure is described in great detail and clarity in Ref. [5lj . 

One can write the differential cross-section as the sum of three terms 

da _ da Iesnm + da pclt da asym 
dPtdy dp\dy dp^ dy dp\dy 

This equation is easy to understand. At low p t , we have the resummed contribution since the latter contributions 
cancel. At high p t we have the perturbative contribution when the resummed and asymptotic cancel. At small p t 
we remove the terms from the perturbative expansion that are asymptotically divergent like 1/pf. This allows for a 
smooth transition between the two regions at all values of pt- However, extracting the divergent pieces can be quite 
difficult analytically as one must express the differential cross-section in terms of p t order by order. For a 2 — > 1 
process, the first order corrections have 2^2 kinematics and this is relatively simple, but becomes more intractable 
for the higher order corrections. 

In this paper, we are interested in the new coefficient functions, and determining where the distributions peak 
for the different colliders, so this treatment will be ignored. However, we will display the perturbative differential 
cross-section to guide the eye on what the transition must look like. 



C. Higgs Resummation 

One of the interesting facets for Higgs production is that there is only a gg initial state for this process in the HQET 
at order af, so the other terms (a qq initial state) are zero explicitly. Without getting too far ahead of our discussion, 
we can see that the Mellin moments of the qq corrections at order strictly vanish on threshold due to the fact that 
there is no qq initial state at lowest order. This makes it possible to work in z-spacc with little additional effort due 
to the presence of the 6(1 — z) terms in the C*gg coefficients, which makes the convolution with the PDFs trivial. 

An additional complication arises from evaluating the parton distribution functions at very low scales during the 
convolution. This is solved by what is known as the 6* t>rescription |5lt l52j . Here the b parameter is replaced by 6* 
that has an infrared cut-off & max so that as b becomes large, 6* — ► 6 max , and the fraction bo/b in Equation 1|29|) never 
leaves the perturbative regime of the parton distribution function. Over the rest of the range b+ ~ b. This can be 
achieved by in the following paramctrization 



K v/l + & 2 /^ax' (37) 

This construction may seem a little artificial, but it allows for the numeric integration of our differential cross-section 
and allows us to use what is known to make reasonable calculations. 

In these calculations, we have set 6 max = (2 GcV) -1 . This is mostly determined by the limits of applicability for 
the PDFs implemented to obtain the hadronic cross-section. There are other non-perturbative correction factors that 
are employed d Ellill, but as no data is available yet we have not included these factors in our analysis. 

It is thus possible to determine the unknown coefficients to a given order in the resummation and preform the 
resummed calculation. The known coefficients for scalar Higgs production will be given later. To leading-log (LL) 
accuracy, only the A^p term is needed. At next-to-leading log (NLL) accuracy one needs the Af \ bP, and C$ 
coefficients. The state of the art currently is NNLL where the B^\ and cf^ coefficients are needed [T^. For a 
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scalar Higgs, these terms have been recently calculated but are missing for a pscudoscalar Higg s. The A g ' term can 
now be determined thanks to the excellent recent work on the three- loop splitting functions |53t l54j . Previously, only 
a numeric estimate was available|55|. 

Let us begin by extracting the process dependent C a b coefficients for the total cross-section. To extract the 
formally divergent pieces of the cross-section, consider the Mellin transform of the hadronic cross-section, u^r(M|). 
The TV— moments in Mcllin space are defined as 



a N {Ml)= f dzz N - 1 a(z,Ml) (38) 
Jo 

The advantage of transforming to Mcllin space is that the limit z — > 1 corresponds to the limit of N — > oo. This 
allows for a systematic way of extracting the divergent terms, which diverge as In (AT) in Mellin space. 

Before continuing, we should comment on which initial state channels contribute. In evaluating the the C a b coeffi- 
cients in the N — ► oo limit we see that only the gg channel has finite contributions, all the other channels have Mellin 
moments are strictly zero on threshold. We could also see that in the HQET there are no qq or qg initial state that 
contribute at this order to the cross-section at p t — 0. Although we have set up our formalism for the sum of several 
channels, we will now consider only the gg initial state channel. 

The Mellin moments of the fixed order corrections allow us to determine the process dependent total cross-section 
C gg ^ coefficients in a simple way. We find the Mcllin moments of the fixed order corrections, A^^, with the package 
harmpol in form|^H|. Some diverge as ln(AT), most tend to zero as N — > oo, and some finite pieces are left over. In 
this way, we can separate the formally divergent pieces from the finite contributions on threshold. With this we can 
identify the pt divergent terms with the ln(iV) divergent terms in the Mcllin moment. We choose to absorb the extra 
powers of 7b into our definition of A^ = Ne lE so that we do not have spurious factors of ^e in our expressions. This 
seems to be appropriate as in the MS scheme the factors of 7# are also absorbed. This being noted, we will continue 
to write our terms as ln(AT) with no factors of je ■ 



1. Next-to-leading-log Differential Cross-section Coefficients 



We are interested in determining the C gg coefficient for the scalar and pseudoscalar Higgs boson. First, let us 
integrate the differential cross-section for the scalar around pt — 0. We will label this contribution 'real' to note that 
this integral is similar to the real emission corrections to the total cross-section. Wc find 



Aa real = 



H ( 47TMr 

° I Ml 



zr(l + e)(l + e)^ 



C A I A) C A i 2 

e 2 e 2 



9t 



(Ml 
{ 4 



(39) 



We have to add the total partonic cross-section virtual correction to this expression to cancel the e poles. The pole 
proportional to the (3q gets renormalized into the coupling like in the total cross-section calculation. Once these two 
expression arc added together we find 



Act 



Oh 

IT 



■^ln 2 



-(f) 



+ C A + -C A + 3C A ( 2 



G 



(40) 



The coefficients can now be read off and agree with the literature 



A^ H = C A , B^ H = -(3o = -(jC A -ln f T R y ^C A + \c A Q 2 . (41) 

As noted earlier, the pseudoscalar Higgs has different virtual corrections from the scalar Higgs boson. This changes 
the Cg^g coefficient for the pseudoscalar to 



C£> A °=C A + IC A { 2 . (42) 

The pseudoscalar coefficient is larger that the scalar coefficient by a factor of 1/12C A . This is a small numeric 
difference, and the NNLL coefficient have not been extracted for the differential distribution, although as we will 

— (2) 

see in the next section wc might expect a larger difference to appear in the differential C gg coefficient based on the 
differences in the C„„ coefficients for the total cross-section. 
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2. Next-to-leading-log Total Cross-section Coefficients 



Exact expressions for the fixed order NLO corrections to scalar and pseudoscalar Higgs production have been in 
the literature for some time. Leaving aside the Sudakov terms (A^ n > and B^ n >), let us examine our expressions for 
the fixed order corrections to the partonic cross-section. We see that the NLO corrections have organized themselves 
in terms of constant pieces proportional to 8(1 — z) from the soft and virtual corrections and additional logarithmic 
corrections. The Mellin moment of the 8(1 — z) is simply 



dzz N ~ 1 S(l-z) = 1. 



(43) 



So it is easy to sec that all the constant terms proportional to 8(1 — z) contribute to the Cgg term. Beyond 
these terms, the Mellin moments of the logarithmic corrections in the limit N — > oo can have finite pieces that also 
contribute. Once the expression for the correction term has been transformed into Mellin space, there are no terms 
proportional to 8(1 — z), but are only constant terms. 

(n) 

In presenting the expressions for the C gg terms, we mix the notation somewhat to allow the reader to see all 
the different contributions. We keep the ln(7V) pieces that are formally divergent, we separate out the terms that 
were initially proportional to the 8(1 — z) for convenience, and we include the terms proportional to ln(M|//x 2 ) for 

(n) 

completeness. We have set \x = fi r = f° r simplicity. The finite pieces compose the C gg coefficients. 
In the case of the scalar and pseudoscalar 



N-gg ~ J 



dz 



-i A (D 

gg- 



H = C A ln 2 (N) 



M 2 

2C A In -4- ln(iV) 



2CUC2 



2CUC2, 



(44) 



A N gg = J lm 



dz 



5 JV-i A (i) 



gg^A 



=C A \n 2 (N) 



2C A ln 



M 2 



\n(N) 



2C A + 2C A ( 2 



(45) 



where the terms in the square brackets are the terms that were proportional to the delta function in the expression 
for the NLO correction in Equations 1|21|) and Ij26(l . We have used the convention of absorbing the extra factors of 

(n) 

that appear in other expressions for the C gg coefficients. 



3. Next-to-next-to-leading-log Total Cross-section Coefficients 

Exact expressions for the fixed order NNLO corrections to scalar and pseudoscalar Higgs production are known. The 
NNLO corrections to inclusive Higgs production have been explicitly calculated and presented for the scalar[3 US 
and the pseudoscalar |35l l3fil | . Although the the color factors have been evaluated in [34l IssT ] . they were found to agree 
perfectly with [3r| once the color factors were evaluated. 

This allowed for the determination of the C gg coefficient for both the scalar and pseudoscalar. The scalar result 
was compared with the literature value [l9l| and was found to be in prefect agreement once the factorization and 
rcnormalization scales were set equal to one another and the spurious factors of je were absorbed. For completeness, 
we present the full expression, leaving the color dependence intact and showing the ln(iV) contributions. 
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(46) 



The term in curly brackets all by itself was the piece proportional to <5( X — z) in the original NNLO correction. 

The NNLO correction used as input for the Mellin moment can be found in Rcf . 36J • One should carefully note that 

there are terms in the coefficient that are proportional to T Rl but not n^T R . These terms come from the higher order 

corrections to the gn effective coupling where only the top quark is included in the derivation of the corrections [TH l30l | . 

( 2) 

One of our novel results is the C gg factor of the pseudoscalar, which is very similar to the scalar case. The difference 
between the two coefficients can be written as 



A (2), A" _ a (2),h 

^N,gg ^N,gg 



C, 



n f T, 
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In 2 (A) 
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In 



'top 



(47) 



(2) 

Finally, as we have the NLO corrections for each of the processes, we can compute the B g ' coefficients for each of 
them. The scalar case matches its value in the literature [l9j and the pseudoscalar result is new. They are 



B^ H + y C2 - jjc) + n f T R C F - n f T R C A Q + ^2 

S(2) ,A° =c *n + n Ca _ 3\ + 1 _ TrC k 
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C F C Al 



(48) 
(49) 



IV. RESULTS AND CONCLUSIONS 



Although we have shown the explicit differences in the scalar and pseudoscalar functions, they are numerically 
quite small. Also, as the differences in the corrections become greater, they are suppressed more in a s , leaving the 
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FIG. 1: The transverse momentum spectrum for the scalar and pseudoscalar Higgs boson at the LHC for \y\ < 2.5. The pt 
distribution peaks at approximately 15 GeV. The resummed curve is the NLL resummation, and the perturbative curve is the 
NLO fixed order calculation. The NLO fixed order calculation diverges in the negative direction at small pt. This piece of the 
differential cross-section is not shown for clarity. These two curves cross at approximately p t = 100 GeV/c and stay very close 
thereafter. 



predominant difference in the resummed cross-section is the same factor of 9/4 that appears in the LO cross-section. 
We are also interested in where the resummed pt distribution peaks at the LHC and Tevatron, so we integrated the 
differential cross-sections numerically for a rapidity \y\ < 2.5. 

We implemented the MRST2002 NLO updated parton distribution functions [HUlI and thc MRST2001 LOparton 
distribution functions j5^] in our analysis as well as the CTEQ 6.1M NLO parton distribution functions [60l l6l| . We 
have taken the renormalization and factorization scales to be identical and set equal to /i 2 = M| + p 2 . It should be 
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FIG. 2: The transverse momentum spectrum for the scalar and pseudoscalar Higgs boson at the Tevatron for \y\ < 2.5. The pt 
distribution peaks at approximately 10 GeV. The resummed curve is the NLL resummation, and the perturbative curve is the 
NLO fixed order calculation. The NLO fixed order calculation diverges in the negative direction at small pt- This piece of the 
differential cross-section is not shown for clarity. These two curves cross at approximately p t = 80 GeV/c and stay very close 
thereafter. 



noted that this choice of scale suppresses the width of the resummation peak as pt grows due to the running of the 
coupling constant becoming smaller as the scale increases, but the effect is not a significant one. The LO cross-section 
"normalization factor" in the resummation formalism also shows p t dependence for the same reason when this scale 
is used. 

Our numerical results for the LHC were created with \fS = 14 TeV and the Tevatron with -/S = 1.96 TeV and 
a Higgs mass of 120 GeV/c 2 . We used an NLO one-loop a s (M z ) = 0.1197 consistent with the MRST2002 NLO 
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FIG. 3: The effects of increasing Higgs mass on the transverse momentum spectrum at the LHC for \y\ < 2.5. The Mh = 
250 GeV/c 2 curve peaks at approximately 23 GeV/c. We can clearly see that the resummed curve peaks at higher pt with 
increasing Higgs mass and that the width of the resummed distribution becomes wider with increasing Higgs mass. 



updated parton distribution functions and a s (Mz) = 0.118 for the CTEQ 6.1M NLO parton distribution functions. 
The pt distributions for the scalar and pseudoscalar Higgs bosons at the LHC are shown in Figure ^ an d Figure [21 
shows the same figures for the Tevatron. It would appear that in both cases the factor of 9/4 difference in the LO 
cross-sections is the dominant difference in the small pt region. The perturbative curves are from the same computer 
code that generated the differential cross-sections in Ref. j^y] ■ 

The average transverse momentum and transverse momentum squared at the LHC for a 120 GeV/c 2 scalar and 
pseudoscalar Higgs boson with p t = — 80 GeV/c are < p t >— 27.5 GeV/c and < pi > x / 2 ~ 32.7 GeV/c with a peak 
value at 15 GeV/c. At the Tevatron the average transverse momentum and transverse momentum squared for the 
scalar and pseudoscalar Higgs boson with p t = — 60 GeV/c are < p t >~ 20 GeV/c and < p\ >*/ 2 ~ 24 GeV/c 
with a peak value at 10 GeV/c. These values should be used with caution as they depend on the PDFs used in the 
analysis. 

We are interested in what happens when a much heavier Higgs boson is considered. This would be the case if one 
were interested in a heavy SM Higgs, the heavy scalar H° in the MSSM, or a heavy pseudoscalar Higgs. We also ran 
our code for a Higgs with a mass of 250 GeV/c 2 . We found a few interesting trends. The peak in the differential 
distribution moved to a higher p t as expected[l(j and the width of the peak became much broader. The width of the 
peak is interesting because it is telling us something about the decay width for the Higgs. The scale of the cross-section 
also dropped considerably as one would expect for a heavier final state particle. As the mass of the Higgs became 
very heavy, it became hard to distinguish a discernable peak in the distribution as it became very wide. As we can 
see in Figure this is a pronounced effect. The resummed curve becomes so broad that it does not cross the fixed 
order differential cross-section until very high transverse momentum. 

There is a great deal of interest in understanding the uncertainties associated with the differential cross-section. 
Although it is common to look at the scale dependence of a total cross-section for scalar Higgs production!^, we 
would like to see how our results are effected by changes in the scale factor [i and the uncertainty in the parton 
distribution functions for the differential cross-section. 

In Figure 0] the upper graph shows the scale dependence of the peak of the distribution when the scale factor is 
varied by a factor of ten. The lower graph in the same figure shows how the entire distribution changes when the 
scale is changed by a factor of two. From this lower graph is it easy to sec that the peak of the distribution has the 
most sensitivity to the scale parameter. We define a prefactor to our renormalization scale to allow it to be varied 
with ease. We define /x 2 = 'X 1 (M% + Pt)- When the scale factor is changed by a factor of ten lower (\ = 0.1) the peak 
increased by a factor of approximately 3.1 and when the scale factor is increased by a factor of ten higher (x = 10) 
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FIG. 4: This is the scale uncertainty in the transverse momentum spectrum for the pseudoscalar Higgs boson. The upper figure 
shows the variation on the differential cross-section at its peak near 15 GeV/c over a scale variation of an order of magnitude. 
The upper figure shows both the MRST and CTEQ parton distribution function. The lower figure shows the variation over the 
whole spectrum when the scale is varied by a factor of two. It is easy to see that the largest scale uncertainty is at the peak 
value. 



the peak of the distribution is lowered by factor of approximately 0.46. When the scale is only varied over a more 
reasonable factor of two, then the peak moves by approximately 25%. The overall scale dependence is very close 
to cts(n) running as expected from the <7q° prefactor in the resummation formalism. Overall, we can see that the 
shape of the distribution is not effected greatly by the change in the scale parameter, only its magnitude is changed 
significantly. 

It is well known that the CTEQ gluon distribution is higher at small x than the MRST sets which can be seen in 
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FIG. 5: The uncertainty due to the parton distribution functions for the CTEQ 6.1M parton distribution functions. The 
resummation is done for the gg initial state only and therefore has the largest PDF uncertainty. A qq initial state would have 
a smaller PDF uncertainty. 

the upper graph in Figure but the effect is quite small. Otherwise, the two distributions are very similar and can 
be considered interchangablc in this analysis. 

In Figure the uncertaity due to the parton distribution functions is shown. At the peak of the distribution, we 
see an uncertainty of approximately 10%. Considering only scale variations of a factor of two would lead us to believe 
that there is still approximately a 35% uncertainty in the differential cross-section at its peak. The uncertainty would 
be slightly lower at other values of the transverse momentum due to the scale /x and larger at higher values of the 
transverse momentum due to the PDF uncertainty. 

In this paper, we have calculated the resummation coefficients for pseudoscalar Higgs boson production for both 
the total cross-section, presenting the B g 2 \ C gg \ and cfg coefficients, and the differential cross-section, presenting 
the C g X g coefficient. We have also shown the effects of increasing the mass of the Higgs boson on the resummed 
differential cross-section and performed an analysis of the uncertainties associated with the renormalization scale and 
the parton distribution functions. 
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APPENDIX A: HARMONIC POLYNOMIALS 

Finding the Mcllin moments of the fixed order total cross-section corrections has been made considerably simpler 
with the HARMPOL package in form[5(|. In order to use this powerful package, it is necessary to express the 
polylogarithmic expressions in terms of harmonic polvlogarithms |r33l l64j . 

Harmonic polylogarithms are defined recursively in three classes for each weight. To make this clear, let us define 
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three functions 

f(-l;x) = T ^—, f(0;x) = -, /(i;i) = -J_, (Al) 
1 + X x 1 — X 

so we can define the weight w = 1 harmonic polylogarithms as 

E(a;x) = [ dx' f(a;x'). (A2) 
Jo 

Thus the first three harmonic polylogarithms can be written explicitly as 

H(-l;x) = ln(l+z), H(0; x) = ln(x), H(l;x) = -ln(l - x). (A3) 

For higher weight harmonic polylogarithms, we need to generalize the notation. The w-dimcnsional vector m w should 
be broken into the first index and the rest of the vector as rfi w = (a,m w -i). This gives us a general expression for 
the rest of the harmonic polylogarithms recursively, 

H(O w ;x) = — j- hi w x, K(m w ;x) — j dx f(a; x ) H(m ra _i; x ). (A4) 

Although it is easy to find the harmonic polylogarithmic expression for the logarithms and dilogarithms, some 
further work is needed for the dilogarithms with quadratic arguments and the trilogarithms that appear in the NNLO 
corrections. The dilogarithms can be simplified in a very straightforward way using well known relationships. To list 
them briefly the most useful expressions are 



-Li 2 (l - x 2 ) = 2[Li 2 (x) + Lia(-x) + ln(x) ln(l - x 2 )] - Q 2 (A5) 

-Li 2 (l - x) = Li 2 (.x) + \n{x) ln(l - x) - C2 (A6) 

Li 2 (a?)= H(0,l;x)=H 2 (x) (A7) 

-Li 2 (-a?) = H(0, -1; x) = H_ 2 (x) (A8) 



Fewer relationships exist for the trilogarithms. It proved to be very challenging to remove three of the trilogarithmic 
expressions simultaneously from the NNLO corrections. The following expressions were derived from the polylogarithm 
literature |6Ri l6fil | and are presented here for future reference (using the notation for the Harmonic polylogarithms of 
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Ref. That allows one to express the NNLO corrections completely in terms of Harmonic polylogarithms, 
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- ln(x) hi (1 - x z ) + - ln 3 (l - x 2 ) - C 2 In 



Li 3 (l - x 2 ) = ln(x) ln 2 (l - x 2 ) + Li 2 (l - x 2 ) ln(l - x 2 ) + ( 3 - 2 



G 



ln- 



l-x 1 



H 2;1 (x)-H_ 2il (x) 



Li 3 (l - x 
ial -- 1 



- ln(x) ln 2 (l - x) + ln(l - x)Li 2 (l - x) + Cs - H 2jl (x), 

Li 3 (-x) +C 2 ln(x) + iln 3 (x), 
6 

Li 3 (x) = H 3 (x), 
-Li 3 (-x) = H_ 3 (x), 

-5i. 2 (l - x) = Li 3 (x) + ln(x)Li 2 (x) + - ln(l - x) ln 2 (x) + Ca, 
-SW-x) =H_ 2 _i(x). 



(A12) 

(A13) 

(A14) 
(A15) 

(A16) 

(A17) 
(A18) 

(A19) 
(A20) 



It was also necessary to linearize all the arguments of the natural logs and to partial fraction the inverse powers of 
1 — x 2 . With these above expressions, it was possible to use harmpol to find the Mcllin moments of the correction 
factors. Although most of these expressions were verified numerically, it should be emphasized that these expressions 
were derived so that they would be valid at x < 1, which is where they would be evaluated on threshold. In some 
regions these expressions would pick up imaginary pieces, but since we are interested in corrections to a partonic 
cross-section our expression must stay real. 

Using these expressions, the NNLO corrections in Refs. £3 HM Eil can be reduced to Harmonic polylogarithms so 
their moments can be easily found. 
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